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Abstract 

We have recently proposed a density functional scheme for calculating the ground-state pair 
density (PD) within the Jastrow wave function PDs of the lowest-order (LO-Jastrow PDs) [M. 
Higuchi and K. Higuchi, Phys. Rev. A 75, 042510 (2007)]. However, there remained an arguable 
problem on the A^-representability of the LO-Jastrow PD. In this paper, the sufficient conditions 
for the A^-representability of the LO-Jastrow PD are derived. These conditions are used as the 
constraints on the correlation function of the Jastrow wave function. A concrete procedure to 
search the suitable correlation function is also presented. 

PACS numbers: 71.15.Mb , 31.15.Ew, 31.25.Eb 
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I. INTRODUCTION 



The pair density (PD) functional theory has been expected to be one of the promising 
schemes beyond the density functional theory.- Recently, we have proposed the PD functional 
theory that yields the best PD within the set of the Jastrow wave function PDs of the 
lowest-order (LO- Jastrow PDs).- The search region for the ground-state PD is substantially 
extended as compared with the previous theory. '^'^ On the other hand, however, there remains 
a significant problem related to the A^-representability of the LO- Jastrow PDs.- 

Let us revisit the problem here. We shall consider an Ai'o-electron system. The Jastrow 
wave function is given by 

"^jixi,- ■ ■,xno) = n f i\^i-'^j\)^ssD{xi,- ■ ■,xno) , (1) 

y^No i<i<j<No 

where Xi denotes the coordinates including the spatial coordinate rj and spin coordinate 
rji, and where Ajy^^, /(|rj — rj|) and ^ssni^i, ■ ■ -, xatq) are the normalization constant, 
correlation function and single Slater determinant (SSD), respectively. The LO- Jastrow PD 
is given by^"^"^ 

7S(rr';rr') = |/ (|r. - r,|)|Sf^ (rr'; rr')^^^„ , (2) 

where j^ssd (^'^'s '^^')Ar=Aro calculated from the SSD. In the preceding paper,- 

we have confirmed that Eq. ([2]) meets four kinds of necessary conditions for the A^- 
represent ability of the PD, and may become "approximately A^-representable" .-'^'^ii^iii'i^ 
However, the possibility of it being A^-representable has not been discussed.- This is an ar- 
guable problem that is concerned with whether the reproduced PD is physically reasonable 
or not. 

The aim of this paper is to discuss the A^-represent ability of Eq. ([2]) and to show the way 
to search the suitable correlation function / (|rj — rj|). The organization of this paper is as 
follows. For the convenience of the subsequent discussions, we first examine the properties 
of the LO- Jastrow PD in Sec. II. The sufficient conditions for the A^-representability of 
Eq. (I2l), which are imposed on the correlation function, will be derived recursively in Sec. 
III. Concrete steps for searching the correlation function that meets these conditions are 
discussed in Sec. IV. Finally, concluding remarks are given in Sec. V. 
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II. PROPERTIES OF THE LO-JASTROW PD 



In this section, we shall discuss the properties of the LO-Jastrow PD. To this aim, the 
properties of PDs that are calculated from SSDs are investigated. The cofactor expansion 
of ^ssoixi, ■ ■ ■, xatq) along the Noth row leads to 

'^ssoixir ■ ■,xno) = ^= { (-l)^°^VA(a;wo)'^'5SD(^ir ■ ^ajjvo-i) 

+ ■■■ (3) 

+ (-l)2^V€(a:iVo)$55D(a;i,- ■ ■,a;^o-i)}, 

where (px^XNo), ^/^(a^Aro), ■■■ and 4>({xno) are the constituent spin orbitals of $5515 (xi, xatq). 
In what follows, suppose that these spin orbitals are given as the solutions of simultaneous 
equations of previous work,- and therefore they are orthonormal to each other. ^^ggj^{xi, ■ ■ 
■,Xno-i) A^o) in Eq. ([3]) denote {Nq — l)-electron SSDs that are defined as the minor 

determinants multiplied by 1/ a/ {Nq — 1)!. The PD that is calculated from Eq. ([3]) is given 
by 

755^(1^'; rr')iv=iVo = iTr ^ Yl ^ssoi^^''^ rr')5v=7Vo-i> (4) 

where 7_^^j^(rr'; rr') 

Ar=Afo-i calculated from ^^ssni-'^i^ ■ ■ '^x^^j-i). 

Likewise, the cofactor expansion of each ^ssoi^i^ ■ x^^^^i) yields (A^^q— 2)-electron SSDs, 
which are denoted by ^^ssni^^^' ' 'iXno-2) (1< j<-A'o — 1). Then, each 755/5 (^'^'s ^^')Ar=Afo-i 
is given by ^ -ifg^{YY'] vv%^^^^_^/{Nq - 3), where -ffg^{vv'] rr')^^;vo-2 is the PD calcu- 

lated from ^'ssoi.^^^' ' Thus, 7^^^(rr'; rr')7v=Afo can be expressed by the sum of 

755^(rr'; rr')w"=Aro-2- the repetition of this procedure, we arrive at 

No No-l 4 3 

7f^(rr'; rr')^=^o = -^^--^ E E " " " E E ^"sU^^''^ ^rO^^lf, (5) 

^ ° 1=1 j = l p=l q=l 

where -ffsoi^^''^ rr')^'^f (l<^<A/'o,l<J<A^o-l,---,l<P<4, l<g<3)are PDs calculated 
from the two-electron SSDs. These two-electron SSDs, which are denoted by ^^^^''(xi, X2), 
are obtained by the above-mentioned successive cofactor expansions. 
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Multiplying both sides of Eq. (j5]) by | / (|r — r'|)| , we finally get 

A^o A'o-l 4 3 



^ ^ ^ i=l j = l p=l q=l 



This relation is the starting point to examine the A^-representability of the LO-Jastrow PD. 

III. SUFFICIENT CONDITIONS FOR THE iV-REPRESENTABLITY OF THE 
LO-JASTROW PD 



We shall start with considering the A^-representability of \ f (|r — r'\)f -y^g^^rr'; rr';^^2 
that appears in the right-hand side of Eq. ([6]). Suppose that the two-electron wave function 
that yields \ f (|r — r'|)|^7^^J^(rr'; rr'Y^^^'^ is given by the following Jastrow wave function; 



--f{\v,-v,\)^%l\x,,x,), 



(7) 



where ^'^ is the normalization constant. From Eq. ([7]), the PD is calculated as 
1/ (|r — r'DI^ 7^^^^(rr'; rr')^'^2''/^2' Therefore, we get A^^""^'^ = 1 as the sufficient condi- 
tion for the A^-representability of |/ (|r — r'|)|^ l^ssni^''^''^ ^^')%=2^- Hereafter, we assume that 
such conditions hold for all values of ■ ■ ■,p and q, i.e., 

j^tr--pq _ ^ g^jj values of i, j, ■ ■ ■,p and q. (8) 

Note that the normalization constant 742^"'^*^ is determined by both /(|ri— r2|) and 
^5S^''(^i) ^2)-^'^ Therefore, if they are given, we can calculate A2 "^'^, and check whether 
the conditions Eq. ([8]) are satisfied or not. As mentioned later, the conditions Eq. ([8]) are 
the parts of sufficient conditions for the A^-represent ability of Eq. ([2]). 
Under the conditions Eq. ([H]), we have 

1/ (|r - r'DI^ ^Psnirr'; rr')^:^ = {^11'^] l^'\rr'; rr%=, l^Zl') , (9) 

where 7'^^)(rr'; rr')Ar=„ denotes the PD operator for an n-electron system. Substituting Eq. 
into Eq. and rearranging, we get 



i=l 



- No~ir , No-2r 5 4 r 3 

-J—y ^— y ...-y -y -yU 

j=i [_ k=i L 0=1 L p=i L 9=1 



ijk-opq 
N=2 



7®(rr';rr')^=2 



ijk---opq 
N=2 



(10) 



It should be noticed that the right-hand side of Eq. (fTOj) has a characteristic form. Con- 
cerning the iV-representabihty of this form, the following theorem holds: 
Theorem. If there exists the set of functions {aa{xn+i)} (1 < a < n + 1) that satisfy the 
conditions; 

a*(x„+i)aa'(x„+i)dx„+i = (11) 

n + 1 

n+l n+1 

a aa{Xn+l)'^%=n{Xli ' ' -^^n) = (- 1)"^ ^ a„ (x„+i) , ■ ■ -^Xn), (12) 
a=l a=l 

then the following equations hold: 

. n+l 

7 E (^^-1 ^^'^(^^'' ^^V=n \n=n) = {^N=n+l\l^'\rv'; rr')^=.+l \<fN=n^l) (13) 

n — 1 ^-^ 
with 

n+l 

^Ar=„+l(Xi, ■ ■ Xn+l) = ^ aa(Xn+l)^^=„(a;i, ■ ■ -, Xn). (14) 

a=l 

Here ^^^„(xi, • ■ ■, x„) (1 < a < n + 1) denote the n-electron wave functions, and o" is a 
permutation operator upon the electron coordinates, and a is the number of interchanges in 
a. 



Proof. The left-hand side of Eq. (1131) seems to be related to the average of ^^(^rr'; rr 
with respect to a density matrix for a mixed state. Indeed, if the density matrix for the 
mixed state is given by 

n+l ^ 

T.^l\'^N=n){n=nl (15) 



lN=n 



n+l 

+ 

then the average of ^^^(rr'; rr')Ar=„ is calculated as 



a=l 

Ml 



, n+l 
' t ~r J- 

a=l 

By using Eq. ( fT6|) . the left-hand side of Eq. ( fT3l) is rewritten as {{n + l)/{n — 
l)}Tr[p„7(2)(rr';rr')^=„]. 

On the other hand, it is expected that the average Tr[p„7^^)(rr'; rr')Ar=„] may be given as 
the expectation value of 7'^^)(rr'; rr')jv=n with respect to a pure state for the whole system 
that includes the n-electron system as a subsystem.— We shall take an {n + l)-electron 
system as the whole system, and suppose that the wave function for such the (n + l)-electron 
system is given by Eq. (fT4l) with Eq. (fTTl) . Then we indeed get 

(v^.+il 7® (rr'; rr')^=n |^n+i) = TT[M^{rr'; rr')^=„]. (17) 



Furthermore, if '^N=n+i{xi,- ■ -jXn+i) is antisymmetric, i.e., if Eq. (fT2|) holds, then the 
expectation value of j^{rr'] rr')iv=n+i with respect to \E'Ar=„_|_i(a;i, ■ ■ ■, Xn+i) is given by 



^„+i| 7(2)(rr'; rr')iv=„+i |^n+i) = (^„,+i| 7^'^(rr'; rr')^=„ • (18) 



n + 1 



Consequently, Eq. f[T3|) immediately follows from Eqs. f[T6l) . flTTI) and f[T8|) . This means 
that Eq. ( fT3ll holds under the conditions that there exists the set of functions {aa{xn+i)} 
(l<a<n + l) that satisfy Eqs. ^ and Q.E.D. 



First, we apply the above theorem to the term y S ( ^jv=2°^^ 

q=l \ 

that appears in Eq. (ITOl) . According to the theorem, this term can be rewritten as 



jijk---opq 
^ N=2 



1 ^ 

/ ^,ijk---opq 

9=1 



7(2) (rr'; rr')Ar=2 



.j,ijk---opq\ 
^ N=2 I 



^,ijk---op 
^ N=3 



7^^-'(rr'; rr') 



N=3 



^.ijk---op 
^ N=3 



(19) 



with 



^f=3°'(^i,a:2,X3) = 5^af-n^3)* 



N=2 V**"!' -^2^) 



(20) 



9=1 



if there exists the set of functions {a^''^ '°''^(x3)} (l<q'<3) that satisfy the following condi- 
tions: 

' ' (21) 



xf ■■■°P(x3)*af -"^^ 



9=1 



ijk---opj 



{x^W^%°^'^{Xi,X2) 



-If J2 af ■••°^(^3)^f=2°"'(^i, ^2). (22) 



q=l 



In addition to Eq. ([8]), the existence conditions for {a^y^''"'°^{xs)} (l<g<3), i.e., Eqs. ([2T 



and ( !22l) . are also the parts of sufficient conditions for the A^-representability of Eq. (12]). We 
assume that the set of functions {a*-''^'°^(x3)} (l<g<3) is obtained. Substitution of Eq. 
(120]) into Eq. ^ leads to 



l/(|r-r'|)|Si^(rr';rrO^=^o 



iVo 



iVo-ir ^ iVo-2r 51- 4 
,-2^ 3^ 2^\ 

0=1 L p=l 



j=l L k=l 



N=3 



^^^\rr'; rr')N=3 



ijk---op 

N=3 



.(23) 



Next, we apply the theorem to the term | ^ \^7v=3°^ 

p=i ^ 



7(2)(rr';rr')Ar=3 



^^^'=3°'' ) that 



appears in Eq. fl23]) . Similarly to the above, we obtain 



,ijk---op 

N=3 



p=l 



7('Hrr';rr')^=3 



,T,ijk---op 
^ N=3 



,T,ijk---o 
^ N=A 



7(')(rr'; rr')^=4 



,T,ijk-o 
^ N=A 



(24) 



with 

4 

nt:{x,,-- ■,x,) = J]af ■■■°(x4)vl/£°^(xi, X2, xs), (25) 
p=i 

if there exists the set of functions {ap^"'°(a;4)} (1 <P<4) that satisfy the following conditions: 

J af-''{x,raf-"{x,)dx =^ V, (26) 

4 4 

aY,af-%x,)^f^,"%x„X2, X,) = (-1)'^ ^ af •••°(x4)vl/j^t3°'(^i, ^2, X3). (27) 

p=i p=i 

The existence conditions for {a'^J^"'°{x4)} (l<p<4), i.e., Eqs. fl26l) and fl27j) . are added to 
the set of sufficient conditions for the A^-represent ability of Eq. ([2]). At this stage, Eq. ([H]) 
and the existence conditions for {ag^''"'°^{x3)^ (l<q'<3) and {a^J''"'°{x4)} (l<p<4) belong 
to the set of sufficient conditions. 

Thus the theorem is applied repeatedly, so that further conditions are added to the set 
of sufficient conditions. Continuing until the conditions for {ai{xNo)} (1 < < ^o), we 
eventually obtain sufficient conditions for the A^-represent ability of Eq. ([2]). 



IV. CONCRETE STEPS FOR CONSTRUCTING THE AT-REPRESENTABLE LO- 
JASTROW PD 

In the preceding section, the sufficient conditions for the A^-representability of the LO- 
Jastrow PD are derived. In this section, we consider the concrete steps for searching the 
correlation function that meets these conditions or checking its existence. 

1. First, we give a trial form of the correlation function. Using this, simultaneous equa- 
tions for the A'o-electron system are solved in a self-consistent way.- 

2. Let us consider the SSD that consists of the resultant spin orbitals for the simultaneous 
equations. The SSD can generally be expanded using the cofactor. The SSD for the 
A'o-electron system is expanded along the Noth row, then we get the A"o number of SSDs 
for the (A^o ~ l)-electron system, i.e., ^ssd i^i, ■ ■ ■, x^^^i) (1 < i < A'o). Successively, 
each of the SSDs for the (A^o ~ l)-electron system is expanded along the (A^o ~ ^)th 
row, and then the (A^o — 1) number of SSDs for the (A^o 2)-electron system, i.e., 
^^SSD (xi, ■ ■ ■,xno-2) (1 < j < Nq — 1), can be obtained for each i. After that and 
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later, the cofactor expansions are likewise repeated, and we finally arrive at the SSDs 
for the two-electron system, i.e.,^^^^^"^^ (xi, X2). The number of two-electron SSDs 
thus obtained is Nq\/2, since j, k, ■ ■ ■■, o, p and q are integers such that 1 <i <Nq, 
l<J<No-l, l<k<No-2,- ■ ■■, l<o<5, l<p<4 and l<g<3, respectively. 

3. The two-electron SSD ^^ssd °^'^ i^T-^ ^2) and the correlation function determine the 
normalization constant A^^''"""^'^:^^ We check whether all of A^J'^""°^'^ are unity or not. 
If no, we return to the step [T] and change the form of the correlation function. This 
process proceeds until all of A^^''"""^'^ become unity. Suppose that such the correlation 
function is found, we get two-electron antisymmetric wave functions ^^=2 °^^(^i' ^2) 
for all cases of i, j, k, ■ ■ ■■, o, p and q through Eq. ([7]). 

4. By means of these ^l/jv=2 "^"^ (^i' ^2); three-electron wave functions are defined as Eq. 
(!20|) . We will check whether there exists the set of |a^''^ " °P(x3), q = 1, 2, 3} that are 
satisfied with Eqs. (1211) and fl22l) . Note that the check has to be performed for all cases 
of i, j, k, o and p. If suitable {a^^^""°P{x3) } cannot be found, we again return to the 
step[I]and modify the correlation function. Suppose the suitable {a^-'^' '"^(xs) } is found 
in this process, then three-electron antisymmetric wave functions ^I//^=3 °^ X2,X3) 
for all cases of i, j, k, ■ ■ ■■, o and p can be constructed from Eq. (!20l) . 

5. We successively proceed the case of four-electron wave functions that are defined 
as Eq. ( l25i) . In a similar way to the step HJ we check whether the set of 
\^a^^'^ ""°{x4), p = 1, 2, 3, 4} meets the conditions of Eqs. fl26l) and fl27j) . If no, we 
restart from the step [1] with the modified correlation function. Suppose that the 
suitable {ap^^ " °(a;4)} are found for all cases of i, j, k, ■ ■ ■■ and o, the four-electron 
antisymmetric wave functions ^E'^=4 °(2^i) ^4) can be obtained from Eq. (125]) . 

6. Likewise, we further proceed the problem constructing the antisymmetric wave func- 
tions for more-electron systems. We search the set of a = 1, 2, ■ ■ ■■, n + 1} 
that are satisfied with Eqs. ( ITTi) and (fT2|) . together with modifying the correlation 
function. If we successfully find the correlation function that meets the conditions 
flTT]) and f[T^ for any n {< Nq — 1), the LO-Jastrow PD of the A^^Q-electron system 
becomes A^-representable. 
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As easily inferred, the above steps are feasible only for the small-electron systems from 
the practical viewpoint. However, it should be noted that the correlation function that 
makes the LO-Jastrow PD A^-representable may, in principle, be found along the above 
steps, though there is a possibility that the suitable correlation function may not exist in 
some system. 

V. CONCLUDING REMARKS 

In this paper, the sufficient conditions for the A^-representability of the LO-Jastrow PD 
are discussed. Using the properties of the LO-Jastrow PD, we derive the sufficient conditions 
that are imposed on the correlation function of the Jastrow wave function. As shown in Sec. 
IV, additional steps to search the suitable correlation function, which satisfies the sufficient 
conditions, are attached to the computational scheme proposed previously.^ Although the 
number of steps rapidly increases with that of electrons, the concrete steps that are presented 
in Sec. IV are feasible for a small-electron system. Of course, there is a possibility that the 
suitable correlation function cannot be found out. In this mentioned in the previous 

paper, LO-Jastrow PDs are approximately A/^-representable in a sense that they satisfy four 
kinds of necessary conditions.->^'^i^»^'ii»i^ 
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